Generalized matrix Ansatz in the multispecies 
exclusion process — partially asymmetric case 



Chikashi Arita 1 , Arvind Ayyer 2 , 
Kirone Mallick 1 , Sylvain Prolhac 3,4 

1 Institut de Physique Theorique CEA, F-91191 Gif-sur-Yvette, France 

2 University of California Davis, One Shields Avenue, Davis 95616 USA 

3 Zentrum Mathematik, Technische Universitat Miinchen, Germany 

4 Department of Physics of Complex Systems, Weizmann Institute of Science, 
Israel 

E-mail: chikashi.aritaacea.fr, ayyer9math.ucdavis.edu, 
kirone .mallickOcea. fr , prolhacSma. turn. de 

Abstract. We investigate one of the simplest multispecies generalization of 
the asymmetric simple exclusion process on a ring. This process has a rich 
combinatorial spectral structure and a matrix product form for the stationary 
state. In the totally asymmetric case operators that conjugate the dynamics 
of systems with different numbers of species were obtained by the authors and 
reported recently [4] . The existence of such nontrivial operators was reformulated 
as a representation problem for a specific quadratic algebra (generalized matrix 
Ansatz). In the present work, we construct the family of representations explicitly 
for the partially asymmetric case. This solution cannot be obtained by a simple 
deformation of the totally asymmetric case. 
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1. Introduction 

The fundamental quest of statistical mechanics is to derive macroscopic behavior from 
microscopic laws. In this respect, the asymmetric simple exclusion process (ASEP) 
has been playing a central role for the last two decades |3S] . The ASEP is one 
of the few models for which the hydrodynamic limit has been proved mathematically, 
using elaborate large deviation techniques [20j [34]. In noncquilibrium statistical 
mechanics, the exclusion process, being one of the simplest non-trivial interacting 
particle process that one can imagine, has reached the status of a paradigm [T{?1 I3"3"] . 
The ASEP is a lattice-gas where each particle is a (possibly biased) random walker 
hopping from a site to one of the neighboring locations only if the target site is 
empty. This exclusion constraint mimics a hard-core interaction amongst particles. If 
the hopping rates are not isotropic, a non-vanishing flow of particles is transported 
through the system which is therefore permanently driven out of equilibrium. Despite 
its simplicity, the ASEP displays a deep mathematical structure that has enabled 
many researchers to perform precise analytical studies. These results can be used as 
benchmarks in the ongoing construction of a general theory of nonequilibrium systems 

ma nu da cna eh eh eh eu . 

The standard exclusion process involves particles and holes hopping on a 
one-dimensional lattice, although many variants have been studied [7J 131) . A 
straightforward generalization of the ASEP is to consider the case of multiple species of 
particles with hierarchical dynamical rules [TJ [5] El 122 • I n fact, such a model with 
A different species of particles automatically appears when one couples N standard 
exclusion processes |18) in a natural way. The A-ASEP provides a fundamental 
example of a multicomponent non-equilibrium process; it has highly non-trivial steady 
states, which are not Gibbs states in general and depend on boundary conditions 

GHEES GI2]. 

For the totally asymmetric case (7V-TASEP), the stationary state was constructed 
combinatorially by Ferrari and Martin |15j . This construction was restated as a matrix 
product Ansatz in [14j and was generalized in [24] to the partially asymmetric case 
(A-PASEP). 

In a recent study [I] of the A-TASEP, the matrix product Ansatz was generalized 
in order to construct a conjugation operator that embeds the [N — 1)-TASEP in the 
A-TASEP. By considering the whole family of A-TASEP processes, with varying 
A, a network of mappings can be constructed (corresponding to an underlying 
partially ordered set -poset- structure). It was shown that the Ferrari and Martin 
construction was a special case of a more general algorithm, corresponding to a 
generalized matrix Ansatz that allows one to lift information from a system containing 
less species of particles to a system containing more species, by recursively splitting 
identical classes of particles into different species. Moreover, the information that is 
obtained is not restricted to the steady states but also affects subsets of the spectrum. 
However, the results of [I] were only valid for the iV-TASEP: the general Ferrari 
and Martin algorithm cannot directly be applied to the iV-ASEP model. Besides, the 
representations of the generalized matrix Ansatz that were obtained for the A-TASEP 
could not be deformed to the A-ASEP in general. The purpose of the present work 
is to fill this gap: we shall explain how to construct all the conjugation operators by 
providing an explicit representation for the quadratic algebras involved in the .ZV-ASEP 
poset structure. We shall also explain how the A-ASEP relates to the Perk-Schultz 
model and investigate the relationship between the conjugation operators and the 
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Perk-Schultz transfer matrix as well. 

The outline of the present work is as follows. In Section [2j we define the model, 
discuss its basic features and briefly review the spectral inclusion properties and the 
generalized matrix Ansatz. The Ansatz starts with a local relation which gives a 
quadratic algebra. The representation for this algebra in the totally asymmetric case 
[4], however, cannot be extended directly to the partially asymmetric case (the N- 
PASEP). We will present a family of representations for the PASEP in Section [3] 
In Section |4j we explain the relevance of our results to the Perk-Schultz model. We 
present the conclusion of the present work in Section [5] 



2. Known results about the A-ASEP model 



2.1. Definition of the model 

Consider a ring with L sites, where a variable (local state) fej <G {1, . . . , N + 1} is 
assigned to each site i £ 1^. Nearest neighbor pairs of local states JK are interchanged 
JK — > KJ with the transition rate 

with 9(0) = 0. Without loss of generality, we set < q < p = 1. We say that the 
site i is occupied by a Jth class particle if kt = J < N. We regard the site i as being 
empty if ki = N + 1. When q = 0, the model is totally asymmetric and is called 
the iV-TASEP; for < q < 1, the model is partially asymmetric and is called the 
N- PASEP. The special case q = 1, corresponds to the symmetric simple exclusion 
process (SSEP). 

The dynamics of the A-ASEP can be encoded in a master equation. Let 
{|1), . . . , \N + 1)} be the basis of the single-site space C N+1 and \ki . . . k^) be the 
tensor product |fci)(g> • • • ®|fci) e (£. N+1 )® L . In terms of the probability vector 

\P(t))= ]T P(k 1 ---k L ;t)\k 1 ---k L ), (2) 

l<fci<7V+l 

the master equation that governs the system is expressed as 

||P(*)) = MW|P(i)). (3) 
The linear operator AfW has the form 

MW = £ (iffl . (4) 

— \ / ^,^+l 



where the local operators [M^2) act on the ith and the (i + l)th components 

V / i.i+l 

of the tensor product and are given by 

N+l 

E (-®(J-K)\JK)(JK\ + e(J-K)\KJ)(JK\). (5) 

J,K=1 

The Markov matrix has zero as an eigenvalue and the corresponding eigenvector is 
called the stationary state. For a given number m, e N of particles of type i (with 
1 < i < N + 1), the stationary state is unique. All the other eigenvalues have strictly 
negative real parts, which characterize the relaxation to the stationary states. The 
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Markov matrix defines an integrable model that can be solved by means of the nested 
algebraic Bethe Ansatz [T][5]. Besides, in [5], the spectral structure of the TV- ASEP 
Markov matrix was investigated; one simple result is that the spectrum of the Markov 
matrix of N species of particles contains the spectra of systems with N' < N species. 



2.2. Matrix Ansatz for the stationary state 



Although the A-ASEP is integrable by Bethe Ansatz, the explicit calculation of 
eigenvectors is not easy even for the stationary states. An alternative technique for 
studying the stationary state is the matrix product Ansatz (which we call simply 
matrix Ansatz), first used in |13j . This trick has grown into very powerful method 
to analyze one-dimensional systems out of equilibrium [7] [TT] . The matrix Ansatz for 
the stationary state of the A- ASEP was constructed in P31 [21] . The basic idea is to 
write the stationary weight for a configuration k\- ■ ■ k^ as a trace of a matrix product 



(6) 



where the operator corresponds to the fc-th species and Z is a normalization 
factor. The operators X^(k = 1, . . . , A+l) must obey specific relations in order for the 
expression ^ to represent the stationary state. We emphasize that the representation 
space of these operators is not the physical space but an abstract vector space which 
is usually infinite dimensional. It was shown in [14] [24] that the X^s can be chosen 
as sums of tensor products of S q ,e q , A q and 1, which satisfy the relations 

S q e q - qe q S g = (1 - q)l, S q A q = qA5 q , A q e q = qe q A q . (7) 

This algebra is related to the quantum harmonic oscillator and the quantum group 
[3] [22l [2~Zl [30]. An explicit representation of this quadratic algebra is given, for 
example, by the following infinite dimensional matrices [13) : 

/ Cl \ / ^ (I \ 



Cl 

c 2 




V 



•7 



Cl 





C2 



A q = 



•7 



(8) 



and 1 = Ai with q = y/l — q ■ 



2.3. Poset structure and spectral inclusion 

Since the number of each class of particles is conserved, the total Markov matrix M^ N ' 
Q splits into blocks as 

M W = M m , with M m e EndV m , where V m = C\h ■ ■ ■ k L ) .(9) 
m #{<|ki=3'}="*i 
We labeled each diagonal block M m and the corresponding vector space (sector) V m 
by m — (mi, . . . , mjv+i)- We will call the label itself a sector, and in particular, a 
basic sector corresponds to the case m, > for all i. A useful alternative labeling of 
the basic sectors [5] is obtained as follows: let Sj be the total number of particles of 
classes k <j, 

sj = mi + m% H h nij. (10) 

One has mj — Sj — Sj-i > (sq = 0), and thus each basic sector can be labeled by the 
set s = {si, . . . , sn} C {1, 2, . . . , £-1} = fl with < s± < s 2 < ■ ■ ■ < s N < L. The set 
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s is an element of S, the power set (the set of all subsets) of ft. In the following, we shall 
use both labels equivalently: for instance, the invariant vector spaces (respectively the 
Markov matrices acting on them) will be denoted either by V m or V 5 (respectively M m 
or Ms). The set S is equipped with a natural poset (partially ordered set) structure 
with respect to the inclusion C, which is encoded in the Hasse diagram. In our case 
it is simply the L — 1 dimensional hypercube where each vertex of the hypercube 
corresponds to a sector, and each edge corresponds to an arrow t —¥ s meaning that 
t C s and #s = #t + 1. (See figure [Tj) 



(1,1,1,1 )++{l,2,3}] 




(4)^0 



Figure 1. The Hasse diagram for L = 4. Each basic sector is labelled in two 
ways. 



The spectral properties of the Markov matrix on the Hasse diagram were 
investigated in j5]. Given two sectors s = {si < ■ ■ ■ < sjv} and t = s \ {s ni , . . . , s nu } 
connected by a finite sequence of arrows in the Hasse diagram, one could define an 
identification operator from s to t as follows 

(p u : \k 1 ---k L )£V s i-> \k[ ■ ■ ■ k' L ) £ Vi with x' = x - #{i\ui < x} . (11) 

Using tfis, a conjugation relation could be proved: (ft s M s = M t tp ls . This relation 
implies that by applying tp is , an eigenvector \E) S of the sector s with an eigenvalue E 
is either projected to an eigenvector (p is \E) s in sector t or is killed out. The surjectivity 
of ip is leads to the spectral inclusion Spec(M 6 ) D Spec(M t ) . 

2-4- Generalized matrix Ansatz and the hat-algebra 

It should be noted that the action of tpt s loses information by projecting a larger sector 
s into a smaller one t. What is desirable is a conjugation operator ( conjugation matrix) 
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ipsi : Vi — > V s satisfying the opposite conjugation relation 

M s i) st = i) st Mi, (12) 
which can be expressed by the following commutative diagram: 

T/ M * V 

1>.t] fv* (13) 

v t v t 

The form for ip st is nontrivial in general, and its action helps us to construct 
eigenvectors (including stationary states) from lower sectors. 

In [1] a method to construct the conjugation matrix -0 st was introduced by 
generalizing the matrix Ansatz for the stationary state. The basic idea is to write 
each element of ip$t m a form 

0'i • ■■3L\^si\kx ■••&£>= Tr(oj lfcl • • • a jLkL ) (14) 

with matrices a JK (1 < J < N + 1, 1 < K < N). 

Consider for a nearest-neighbor pair of sectors s = {si < • • • < sjy} of TV-species 
and t = s \ {s n } of (N — l)-species. The sector t is obtained from the sector s by 
merging the particles of nth and (n + l)st classes into the same class. It will be 
convenient to write these operators as elements of an operator-valued matrix a such 
that 

a JK = (J\a\K) . (15) 



Following [4], one can prove that the matrix - 0s,s\{s„} defined by equation (14) is a 
conjugation operator if there exists another operator- valued matrix g^ N ' n ' of size N + l 
by N that satisfies the following key relation (hat relation): 

M ioc ( a ® a ) " (a^-Mi^T 1 * = a®a - a®a . (16) 
The elements of this equation give a set of relations, which we call the hat algebra. We 



emphasize that the hat relation ( 16 ) connects the local Markov matrix of the 7V-ASEP 
and that of the (N — 1)-ASEP, and is of a different nature from the relations used in 
p"7l [26] for the usual matrix Ansatz for the stationary state. 

Therefore, the problem of constructing a conjugation operator has been reduced 
to finding realizations of the hat algebra. 



2.5. The N-TASEP 



In the TASEP case, a family of solutions (a, a) = ( a ^ n ) , a^'") ) (1 < n < N) to the 



hat relation (16) was successfully constructed in g], where n is the class being split. 
The operator-valued matrices a' iv ' n ' and 

£(N,n) are 

related by multiplication with a 

diagonal matrix: 

= diag(l, ...,1,0,..., 0)a( w '"> . (17) 



JV+l- 



Thcn the hat algebra ( 16 ) becomes a closed quadratic algebra generated by the 
elements of aS N,n K An explicit representation of the hat algebra is constructed by 
(N — l)-fold tensor products of the fundamental quadratic algebra generated by the 
infinite dimensional matrices 6, e and A that satisfy the relations [T3] 

Se=l, 5A = 0, Ae = 0. (18) 
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We note that they are obtained by taking q = in Q. 

The expressions of the operators a^f^f 1 ^ are given in the following table with 

<y®l® (_1) (gie = 1 and 6(g>t® x ®e = for x < -2. 





1 • • • n-l 


n 


n+l ••• TV 


1 

71—1 


A^-VfgiS®!®^--'^ 


J®1®( JV - / - 1 ) 




n 





•j(8>(JV-n) 




71 + 1 









n + 2 
N+l 




]!®(J-3)^ 
e( ^(JV-.7+l) 


^(K-2) 0fel »(J-K-2) 



3. Generalized matrix Ansatz for the iV-PASEP 



In this section, we explain how to construct representations of the hat algebra and the 
corresponding conjugation operator for the iV-PASEP case. It turns out that finding 
a representation for the hat algebra of the iV-PASEP is not a simple deformation of 
the TV- TASEP case. It will require a construction more involved and rather different 
from that used for the iV-TASEP. 



3.1. The hat relations for PASEP 



For any given value of N, the TASEP solution ( 19 ) can be generalized to the PASEP 
in the two special cases n = 1 and n = N. The case n = 1 corresponds to splitting the 
first-class particles in two sub-classes. The dual case n = N corresponds to splitting 
the holes (labelled as AT's) in the (JV — 1)-PASEP model into iVth-class particles 
and holes (now labelled as (N + l)'s) in the TV- PASEP model. The generalization is 
obtained by the replacement {6, e, A} — > {S q , e q , A q } (|8| in each element of aS N ' n ^ and 



a simple modification of |17|) 

a^ 1 ) =diag(l, <?,..., g)^' 1 ), (20) 

S<^=diag(l l ...,l I «)«W^ J (21) 

so that the hat relation ( 16 1 is satisfied. The hat algebra of the TASEP case is 



deformed by this modification. 

For general n (with 1 < n < N), equations (20) and (21 1 naturally lead us to 
assume the following relation between 3^ N ' n ' and a™"': 



i(N,n) = d (N,n) a (N,n) = di&g(1) 



N+l- 



(22) 



Inserting this equation in the hat relation (16), one finds a closed quadratic algebra 
for the operators qjk- 
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The quadratic hat algebra for the A-PASEP is defined by a set of relations 
between its generators, summarized a^|] 





l:K <K' 


II : K > K' 


A : J < n < J' 


ajK'dJ'K — qdj'KdjK' 

= (1 - q)a JK aj> K ' 


qaj'KdjK' = qajK'aj'K 


„ J < J' < n or 
n < J < J' 


ajK'dj'K — qaj>KajK' 


qaj'najK' - qajK'dj'K 
= (1 - q)a JK aj'K' 


C : J = J' 


ajRajK' = a.iK'ajK 


qa.jKa.jK' = qajK'ajK 


D : J > n > J' 


a J'KajK' — a JK' a J'K 


aj'KajK' — qajK'aj'K 
= (1 - q)a JK aj'K' 


J > J' > n or 
n> J > J' 


ajK'dj'K — a.FKajK' 
= (1 - q)a JK aj> K > 


aj'KdjK' 
= qajK'aj'K 



(23) 



Our aim is to construct an explicit representation for this iV-PASEP algebra. 
At first thought, one would expect that the replacement {5, e 7 A} — > {S q ,e q ,A q } 



in the TASEP solution ( 19 ) provides us with a representation of the hat algebra 
(23 1 for any value of n. Unfortunately, this is true only for n — 1 and n = N. 
1 



For general n 



(1 < n 



(N,n) 



"JK 



TASEP 



V b JK 



< N), we found no perturbative representation a JK — 

starting from the TASEP 



+ <fc 



2 (N,n) 



JK 



for the algebra (23 



representation (19). The heuristic reason behind thisfact can be summarized as 
follows. If one considers the fundamental TASEP algebra T generated by <5, e and 
A then the tensor products <5<X><5, e®e and A<&A also generate the same algebra i.e. 
they satisfy the same quadratic relations. Therefore there exists a simple coproduct 
operation from T to T®T that preserves the algebraic structure (181 [5j. However, 
one can easily verify that 5 q (£>5 q , e q ®e q and A q ®A q do not satisfy the fundamental 
PASEP relations Q. It seems that no coproduct exists for PASEP algebra that would 
allow us to build natural tensor representations. This is the mathematical obstruction 
that prevents us from constructing the N- PASEP hat algebra by deforming the known 
iV-TASEP hat algebra. 



To summarize, the solution for the N- TASEP ( 19 1 is of no help in general to find 



representations of the N- PASEP algebra ( 23 1 . One needs a different strategy to build 



explicit representations, which we will explain beginning with a specific example in 
the next subsection. 



3.2. The simplest non-trivial example 

In this subsection, we work out the simplest non-trivial example with 1 < n < N: this 
is obtained for the case (N,n) — (3, 2). In particular, we show that the relations (23) 
are not contradictory by giving an explicit representation. Hence they define a bona 
fide algebra. 

A solution to (23) will be constructed by using the regular representation of this 
algebra. We assume that 

an = a 43 = Id, a 2 i = a 33 = 0, (24) 



| The cases A-II and C-II vanish when q = 0. 
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H13CJ12 


— 


O12O13 




Ol30l2 


- 


012013 








ffll2d22 


= 


0O22O12 




012022 


= 





O12 


-< 


a 22 


0113022 


= 


0022013 




Ol 3 022 


= 





Ol3 


-< 


(122 


O23O12 


= 


O12O23 + (1 ~ 


- g)a 2 20i3 


O2 3 0l2 


= 


012023 + 022013 








O13O23 


= 


0O23O13 




Ol 3 023 


= 





ai 3 


-< 


023 


023022 


= 


O22O23 




023022 


= 


022023 








£112031 


= 


0031012 + (1 


- g)a 32 


O12O31 


= 


O32 


012 


-< 


a-.ii 


ai3CJ31 


= 


ga 3 ioi 3 




013031 


= 





013 


-< 


a-M 


fl22tl31 


= 


0O31O22 




022031 


= 





"22 


-< 


a-.ii 


123031 


= 


0031023 




023031 


= 





023 


-< 


031 


012032 


= 


032012 




O12032 


= 


032012 








013032 


= 


qa 3 2 0i 3 




Ol 3 032 


= 





ai 3 


-< 


"32 


022032 


= 


032022 




022032 


= 


032022 








023032 


= 


0O32O23 




023032 


= 





023 


-< 


O32 


032031 


= 


031032 




032031 


= 


031032 








012041 


= 


ga 4 iai2 + (1 


— 0)042 


O12041 


= 


042 


ai2 


-< 


0,41 


O13O4I 


= 


ga 4 iOi 3 + (1 


-q)Id 


013041 


= 


Id 


013 


-< 


041 


022041 


= 


9041022 




022041 


= 





a 2 2 


-< 


0,41 


023041 




Q041023 




023041 







023 


-< 


041 


031041 




9041031 




03 1041 







a-.ii 


-< 


041 


032041 




0041032 




032041 







032 


-< 


041 


012042 




042012 




O12042 




042012 








O13O42 




9042013 + (1 


- g)ai2 


Ol 3 042 




Ol2 


ai 3 


-< 


042 


022042 




O42O22 




022 042 




042 022 








023042 




0O42O23 + (1 


— 5)022 


023042 




022 


023 


-< 


042 


042 031 




031042 + (1 - 


- 0)041032 


042031 




031042 + 041032 








032042 




0042032 




032042 







032 


-< 


042 


O42 O41 




O41O42 




O42 041 




041042 









Table 1. The hat algebra j23| for (N,n) = (3,2). The left and middle columns 
correspond to the PASEP and TASEP cases, respectively. In the right column 
the ordering restricted by the algebra of the TASEP case. 



as in the TASEP solution. We also assume that equation (22) is valid so that 
a( 3 > 2 ) — (l. l ; g, g) a' 3,2 ) . Inserting these assumptions in the hat relation (16 1, one 
obtains 28 relations, shown in the leftmost column of table [l] that have to be satisfied 
by the 8 unknown a^ 's. Let us consider the space of monomials generated by ay's. 
The TASEP algebra (the middle column of table [lj tells us the correct order of the 
unknowns ay in an arbitrary word consisting of these monomials. 

For example, 012 should be located to the right of 022 (022 >- 112) due to the 
relation 012022 = 0. We listed the restriction for the reordering of the generators in 
the right column of table [T] One of the allowed orderings is 

041 >- «3i >- 042 >- a 32 >- a 22 >- a 12 >- a 23 >- a i3 . (25) 
Let us consider the right-action of each generator ay on the monomial (word) 

W = al?al?---a n x f (ny > 0), (26) 
for the TASEP case. We set 

am = t® 7 ®5 . (27) 
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From the relation 013023 = 0, one possibility for 023 is 
o 23 = 1® 6 ®<5®A 

From the first, forth and sixth relations, the action of an is calculated as 

(nia - 0) 



Wa 12 = 



Ml 

,"41 

Ml 



„™22„ «23 ™13 

■ 0-22 0i20 23 a 13 

„"22 „"12 + I„n23„«13 

' "22 "12 "23 u 13 



'41 



"32 „™ 2 2 + l„«23-l„«13 



' a 3 2 22 



23 



a 



13 



(ni2 > 1) 



Thus one can set 



a 12 = 1® 5 ®<5®1^ + l® 4 ®<5®A®e®<5 



n<8>2 



£4. 



In a similar way, we obtain conditions for the other monomials in the order 022, ■ 
and we end up with 

a 22 = 1® 4 ®<5®A®1®A, a 32 = 1® 3 ®<5®1® 2 ®A® 2 , 

042 = !® 5 ®5®l®e + !® 4 ®5®A®e®l + l® 2 ®^®!® 2 ®^® 2 

a 31 = 1®<5®1® 2 ®,4® 4 + <5®A®e®<5®A® 4 + !® 3 ®5®l®e®,4® 2 

a 4 i = l® 7 ®e + l® 2 ®<5®A®l®e®,4® 2 + 5®A®t®A® 5 . 



(28) 



(29) 

(30) 

j 041 1 

(31) 
(32) 
(33) 
(34) 

Accidentally (and fortunately), the q-deformation of <J's, e's and A's in ay's gives 
a solution to the hat relation for general q. Furthermore, the following simplification 
still keeps the hat relation satisfied. Erase the 1st, 2nd, 3rd, 4th and 5th components 
of the tensor products in each a%j. Erase the 3rd term in 042, the 1st and 2nd terms 
in d3i and the 2nd and 3rd terms in 044. Then we get a solution 



,0,2) 



/ 1®1®1 


e®A®A 
\ l®ll®e 



A®e®8 + <5®1®11 

A<E)t<E)A 

t(g)A<g)A 
A®e®l + <5®l®e 



1®1®<5 \ 
1®5®A 


ligil® 11 J 



(35) 



It is straightforward to verify that the cljk's satisfy the algebra (23). Using this 
representation, one can calculate all matrix elements (14) of the conjugation operator 
^{1,2,3}, {1,3} between the sectors s = {1, 2, 3} and t = {1, 3}. For example, one has 

(1324|V { i, 2 ,3},{i,3}|2132) = Tr(4f } 4 3 2 2) ) 
= Tr ((A®e®6 + <5®l®l)(e®A®yl)(l®(5®A)(A®e®l + (5®l®e)) 
= Tr(AeA) Tr(eAde) Tr(<L4 2 ) + Tr(AeS) Tr(eAS) Tr(<L4 2 e) 

+ Tr(SeA) Tr(A6e) Tr(A 2 ) + Tr(SeS) Tr(Ae) Tr(A 2 e) 
_ 1 + g 2 

~ (l-g2)2(l_ g 3)- 



(36) 
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The full elements of matrix "0(1, 2. 3}, {1,3} are given by: 



^{1,2, 3}, {1,3} 



qb 
b 



qb 



cq qd 3 ■ qd 3 q 2 e \ 

qb q 2 d qe ■ qe qd 

a b ■ b c q 2 d q 2 d ■ q 2 e 

qb a ■ qe b qd ■ ■ q 2 d qe 

qb a qd ■ cq ■ qb e qd 

b-ad-c-b-ed- 

b q 2 d qe qb qe qd 



where we have set 
1 

a 



qb qd a e 

qd 3 qb q 2 e a 

q 2 e ■ q 2 d a 

qe ■ qd qb 

qb qe qd a b q 2 d 

cq qb qd qb a qd ■ e 
b d c a e b 

q 2 d b c q 2 e a q 2 d 

cq qd 3 qb qd 3 q 2 e qb ■ a 

q 2 d qe b qe qd 

qd 3 q 2 e ■ qb 
q 2 d ■ q 2 e b 
qe qd 
e ■ qd 

q 2 d qe ■ qe qd 

e d ■ d c 



1 



d b d 
cq qb 
cq qd 3 

b q" d b 

a q d qe 

qe 



qb 
c 
b 



qd. 
b 



qb ■ a 

cq ■ qd 3 

■ q 2 d ■ 

q 2 d qb qe 

cq qb 



b 
a 
qb 



qb 
b 



qb 

b 



l + q 2 



d = 



(1- 3)3(1 - g a)> 
1 



(I -q)(l -J*)*' 
1 



(l-g2)2(l_g3)' 



(37) 



(38) 
(39) 



(1-^)2(1-33)' (l-g)(l- (? 2 )(1 _ (? 3 ) ' 

and replaced by a dot to make reading easier. The ordering of the bases is 
lexicographic: 1234,1243,. . .,4321 for the sector s = {1,2,3}, and 1223,1232,. . .,3221 
for t = {1,3}. One can check explicitly that the conjugation relation 
"0{i,2,3},{i,3}-^{i,3} = -^1,2,3} ^{1,2,3}, {1,3} is satisfied. We remark that all the nonzero 
elements of this example have a singularity at q — 1 of order 3. 



3.3. The general case 



One can construct a family of solutions to ( 16 1 for the general case (N, n) recursively 
by using the case (N — 1, n) if n < N or the case (N — 1, n — 1) if n = N: 



(4,1) (4,2) (4,3) (4,4) 



(40) 



(3,1) (3,2) (3,3) 

\ \ / 
(2,1) (2,2) 
\ / 
(1,1) 

More explicitly, for 1 < n < N — 1, a^ N,n > is defined in terms of a^ -1 '"' by the 
following formula: 
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a JK - 

E 

J<j<n 

t®( E - J ^®5®A®( J 



a 



1). 



a 



(N-l,n) 
JK 



l<j<n 



5^~ 1,n) ®l® (n - ; ' ) ®e®l® (, '- 1 - J) ®«®A®^- :l J (l<J<n,l<Jr<JV- 

(1 < J< n,Jf = N), 
(n + 1 < J < N, 1 < K < N - 1), (41) 
(n + 1 < J < iV, if = iV), 
(J = AT + 1,1 < AT-1), 



(J = AT+1,# = AT), 



where E = E{N,n) = nN — n 2 + n — 1 is the number of the tensor products, and 



31' 



s(-i), 



)(5 = 1. For n = N ', a^ N ' N ^ is defined in terms of a^ N 1,JV 1 ' 



l JK 



= < 





As 



(JV-l.iV-l) 
1 N,K-1 

(N-l.N-l) 
J-l.K-l 



^18(^-2) 

{N-l,N- 
^ a N,K-l 



1) 



(J = if = 1), 

(J = 1,2 < JK" < iV), 

(2 < J< N,K = 1) , 

(2 < J < AT, 2 < X < N), 

(J = N + 1,K = 1), 

(J = N + 1,2 < K < N). 



(42) 



We remark that in this representation the number of tensor products grows 
typically as iV 2 (supposing that n is of order N) whereas in the TASEP case the 
representation (JI9]) involves only (AT — l)-fold tensor products. Since the trace of each 



component of the tensor product gives the factor , each nonzero element of ip has 
a singularity at q = 1 of order E{N, n). 

This recursion has been obtained by investigating the case (AT, n) = (4, 2) (see 
subsection 3.2 1 and by guessing the general (N,n) case. The general formulae (41 1 
and ( 42 ) are proved by verifying that the quadratic relations ( 23 ) are satisfied. This is 
done by induction on the number of species AT. I.e. one can check all the cases in the 
table (23) assuming these are satisfied for A^ — > N — 1. The proof is not particularly 
illuminating, which is in the same spirit as the one given in [24] . 

When g = and n = 1, N, the operators obtained from (41 ) and (42 ) are identical 
to the TASEP solution (19 1. This is not true anymore for n^l,N. In other words, 
in the TASEP case (q = 0), one has two different families of representations for the 
relations (23). We expect, however, that the conjugation matrices constructed by two 
different representations are identical: 

1ps,s\{s n }\q=Q = i>s,s\{s n } (43) 
for 2 < n < N — 1 as well as for n — 1, N. So far this identity has been checked for 
small systems by using Mathematica. 



3.4. Conjugation paths through the Hasse diagram 

We finally consider the case of two PASEP models with respectively A^ and AT' 
species and we suppose that AT and AT' are not consecutive integers. We show that a 
conjugation operator between the A^-PASEP and the A^'-PASEP can be constructed 
by using recursively the results of the previous sections. 
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First, we define a ^-product on operator valued matrices. Consider two operator 
valued matrices A — {A^}^ and B = {Bij}^ and let the symbol * denote the product 

A*B = {Aij®B jk } ik , (44) 

which is bilinear and associative: (A * B) * C = A * (B * C). When A or B is a 
scalar-valued matrix, ★ is just the usual product. The following formula is satisfied: 

(A®B)*(C®T>) = (A*C)®(B*V). (45) 

We now suppose that we have found solutions for the hat algebra for all N £ N: 

MW(aW®aW) - (a^W^M^" 1 ) = a^Sa™ - a^W* . (46) 

Then, using the ^-product one can construct a solution of the general hat-algebra 
defined as follows: 

M Loc } {X®X) = {X®X)M[ N 'J + X®X - X®X . (47) 

Indeed, one can show by induction that this relation is satisfied by the choice 

^/'uMi.^a^ 1 ) (48) 

X= aW*..^a^+ 1 )*a«*a( l - 1 )*...*a( Ar '+ 1 ). (49) 

N>i>N' + l 



Generically there is no linear relation of the type (22 ) relating X and X ; hence, a 
closed algebra cannot be defined by the elements of X alone (except for the TV' = N— 1 
case). Furthermore, we observe that for 1 < n < n' < N, a^ N ' n ' * a( Ar_1 ^™) ^ 
and therefore X depends on the choice of the intermediate values 



of the n's. In other words, one has several solutions to (47 1. We conjecture, however, 
for a iV-species sector s = {si, . . . ,Sj\r} (si < ■•• < s^,l < n < n' < N), the 
"commutation relation" holds up to an overall factor as 

Y\ N , (1 - g s '- s ") 

4>s,s\{s„}A\{s n }, S \{s„,s n ,} = J7n-tl . , „.s ^^\{» n /}^\{» w >},»\{»n.»v}' ( 50 ) 

or equivalently the following diagram commutes up to an overall factor: 

^ S , S \{a„} / / \/> S , S \{s n ,} 

V >\M V *\{° n '} ■ ( 51 ) 

^\{s„} > s\{ 8 „ 1 v}\ / V's\{ ;s „,}, S \{ s „, s „,} 

V*\{s n ,s n ,} 

So far this identity has been checked for small systems by using Mathematica. 



For N' — 1 the relation (47) corresponds to the usual matrix Ansatz for the 
stationary state (which is unique up to an overall constant in each sector), and 
one has several matrix representations for the stationary state according to the path 
— » • • • — > 5 on the Hasse diagram. The stationary state was first obtained in [21] , 
where the path corresponds to — > {si} — > {s\, s 2 } s \ {sn} — > s. 

Remark: In the symmetric exclusion case p = q = 1 (SSEP), the Markov matrix 
obeys the detailed-balance condition, and is a symmetric matrix: M T = M. Thus 
the system in each sector converges to an equilibrium stationary state, where all the 
possible configurations are realized with an equal probability. Since the conjugation 
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matrix ip has singularity (1 — q) E ( N > n ) ; we take the limit q — > 1 after multiplying it 
by (l-q) E (N,n). 

^ )A{S „ } = lim(l - q) E ^^\ {Sn} (52) 

for s = {si, . . . , sjv} and 1 < n < N. This matrix satisfies the conjugation relation 
dl2b: M s tJj s , s \{s n } = - 0s.s\{s„}^5,\{ s >>}- SSEP case, one can easily find a few 

other possibilities for a conjugation matrix satisfying ( 12 1. For example, the transpose 



of the identification operator tp also satisfies Mgip^ = (p^Mi for s D t since MJ- = M s 
and M t T = M t . Another simple solution (one dimensional representation) to the hat 



l JK 



relation p6|) with (|22|) is 

(2 < J <n,l< K < J -1 

or n + 1 < J < N, J < K < N), (53) 
1 (otherwise). 

Then the matrix ip' defined as (Ji • • ■ JlW s ^\K\ ■ ■ ■ Kl) — Ili<i<L a 'j k a ^ so satisfies 
4. The A-ASEP as a multistate vertex model 

Let us consider the two dimensional vertex model, e.g. on t x L lattice. As 
in figure [2j each vertex has a Boltzmann weight W% t j defined by values 

(frij-ii c i-i,j) c i,j) that can be assigned to its four connected edges. The partition 
function is given by 

z = e n (54) 

configuration vertices 

The Perk-Schultz model defines a family of vertex models that are exactly solvable 
[2"5] . As we will shortly review, the A- ASEP can be realized as a special case of 
the Perk-Schultz models for a particular choice of the vertex weights. We will also 
investigate the conjugation relation for this particular restriction of the Perk-Schultz 
model. 

Define a matrix R(X) S End ( (C Ar+1 )® 2 ) as 



R(X) = p (l + AMffl) , (55) 

where p is the permutation matrix: p(\a)®\P)) = |/3)®|a). By regarding each 
element R™(X) — (xy\R(X)\zw) as the Boltzmann weight of each vertex, the N- 
ASEP corresponds to a special case of the Perk-Schultz model with N + 1 states, 
where the non-zero elements are given as 

a a P 

P a 



a 



(56) 



a P P 

flSS(A) = l, R^(X) = XQ(a-P), R^(X) = l-XQ(a-P), 
see [51)1 [3!] . The matrix R(X) satisfies the Yang-Baxter relation 

R bc (v)Rac(v)Rab(X) = R ab (X)R ac (fi)R bc (v) (57) 

with A = i_( p +q}„ +pqiil/ ■ The indices a,b and c specify the spaces on which the R 
matrices act. The monodromy matrix is defined as 

T a (X)=R aL (X)---R al (X) (58) 
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c m c 02 C03 



c i)L 



C £l C 12 C £3 



C £L 



-631- 



L 22 ^23 



c 2£ 



L 12 ^13 



-611- 



C 01 C ()2 



Cox 



Figure 2. The two-dimensional vertex model. The 7V-ASEP corresponds to the 
case where Boltzmann weights are given as equation (|56[>. 



acting on C N+1 ® (£ N+1 ) where the first space is the so-called auxiliary space and 
is denoted by a. It satisfies the "global Yang-Baxter relation" 

%{v)T a {^)R ab {X) = R ab {X)T a {ii)T b {v) (59) 
with A = l-fp+^T+pq^iy • The row-to-row transfer matrix T(A) is defined by the trace 
of the monodromy matrix over the auxiliary space, 

T(A) = Tr Q T a (A). (60) 

The Yang-Baxter relation implies that it constitutes a one-parameter commuting 
family: 

[T(Ai),T(A a )] = 0. (61) 
The Markov matrix can be rewritten in terms of T(A) as 



(62) 

which indeed satisfies \M^ N \T{X)] — 0. The transfer matrix T(A) and the Markov 
matrix can be diagonalized by using the Bethe Ansatz pQ [51 129j . 

The number of each species of particles is invariant under the action of the transfer 
matrix T(A) as well as the Markov matrix (T(X)V S C V s ) and we denote the restriction 
on the sector s by T S (X). We conjecture that the transfer matrix also satisfies the 
conjugation relations with the identification matrix and with the new conjugation 
matrix up to order X L : for sectors s = {si < ■ ■ ■ < s^} and t = s \ {s n } 

^ 5 T 5 (A) - T t (A)^ = q Sn \ L (fii S , (63) 

T s (\)ip st - i> st T t (\) = q s -X L i> si . (64) 
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Remark: In the nested algebraic Bethe Ansatz technique, eigenvectors of M^ N ^ are 
constructed by the action of a product of "B-operators" (elements of T(A)) on the 
vector |1 • • • 1). On the other hand, the generalized matrix Ansatz also enables us 
to construct an eigenvector by the action of a product of tfi's on an eigenvector with 
the same eigenvalue in a lower sector. In particular, the stationary state can be 
written as if) ■ ■ ■ ip\l ■ ■ ■ 1). It seems that our conjugation operator and the B-operator 
play similar roles. However the generalized matrix Ansatz never gives us information 
about eigenvectors with new eigenvalues in each sector. 

5. Conclusion 

We have applied the generalized matrix Ansatz to the multispecies ASEP, where the 
central issue is reduced to finding a representation for the hat algebra. The family of 
representations that we find here is defined by a recursion formula with respect to the 
number of species N. This new solution is not obtained by perturbation of the known 
solution for the TASEP. The SSEP is also a special case where a one-dimensional 
representation exists. We conjecture that the conjugation relation continues to hold 
(modulo A L ) for the Perk-Schultz transfer matrix. 

As we remarked in the last section, our generalized matrix Ansatz does not 
enable us to construct eigenvectors with new eigenvalues in each sector, which is a 
problem to be solved. Another open question is whether there exists another family 
of representations to the hat algebra. (In fact we have found different families of 
representations for the TASEP and SSEP cases.) Applying the technique to open 
boundary conditions with injection and extraction of particles or other driven-diffusive 
processes with rules more general than equation can also be an interesting study 

0S- 



Acknowledgment s 

It is a pleasure to thank C. Banderier, M. Bauer and G. Duchamp for very useful 
discussions on the deformations of quadratic algebras. We also thank S. Mallick for a 
critical reading of the manuscript. C. Arita is a JSPS Fellow for Research Abroad. 

[1] F. C. Alcaraz and R. Z. Bariev, 2000, Exact solution of asymmetric diffusion with N classes of 
particles of arbitrary size and hierarchical order, Braz. J. Phys. 30, 655. 

[2] F. C. Alcaraz and V. Rittenberg, 1993, Reaction-diffusion processes as physical realizations of 
Hecke algebras Phys. Lett. B 314, 377. 

[3] C. Arita, 2012, Remarks on the multi-species exclusion process with reflective boundaries, J. 
Phys. A 45, 155001. 

[4] C. Arita, A. Ayyer, K. Mallick and S. Prolhac, 2011, Recursive structures in the multispecies 

TASEP, J. Phys. A 44, 335004. 
[5] C. Arita, A. Kuniba, K. Sakai and T. Sawabe, 2009, Spectrum of a multi-species asymmetric 

simple exclusion process on a ring, J. Phys. A 42, 345002. 
[6] R. J. Baxter, 2007, Exactly solved models in statistical mechanics, Dover. 

[7] R. A. Blythe and M. R. Evans, 2007, Nonequilibrium steady states of matrix product form: A 

solver's guide, J. Phys. A 40, R333. 
[8] L. Cantini, 2008, Algebraic Bethe Ansatz for the two-species ASEP with different hopping rates, 

J. Phys. A: Math. Theor. 41, 095001. 
[9] V. Chari and A. Pressley, 1995, A guide to Quantum Groups, (Cambridge University Press, 

Cambridge). 



Generalized matrix Ansatz in the multispecies ASEP 17 

[10] B. Derrida, 1998, An exactly soluble non-equilibrium system: the asymmetric simple exclusion 

process, Phys. Rep. 301, 65. 
[11] B. Derrida, 2006, Matrix ansatz and large deviations of the density in exclusion processes, 

Proceedings of the International Congress of Mathematicians Madrid Spain 2006 367-382. 
[12] B. Derrida, 2011, Microscopic versus macroscopic approaches to non-equilibrium systems, 

Proceedings of Statphys 24, J. Stat. Mech.: Theor. Exp. P01030. 
[13] B. Derrida, M. R. Evans, V. Hakim and V. Pasquier, 1993, An exact solution of a ID asymmetric 

exclusion model using a matrix formulation J. Phys. A 26, 1493. 
[14] M. R. Evans, P. A. Ferrari and K. Mallick, 2009, Matrix Representation of the Stationary 

Measure for the Multispecies TASEP, J. Stat. Phys. 135, 217. 
[15] P. A. Ferrari and J. B. Martin, 2007, Stationary distributions of multi-type totally asymmetric 

exclusion processes, Ann. Prob. 35, 807. 
[16] O. Golinelli, K. Mallick, 2006, The asymmetric simple exclusion process : an intcgrable model 

for non-equilibrium statistical mechanics, J. Phys. A: Math. Gen. 39, 12679. 
[17] H. Hinrichsen, S. Sandow and I. Peschel,1996, On matrix product ground states for reaction- 
diffusion models, J. Phys. A: Math. Gen. 29, 2643. 
[18] T. M. Liggett, 1985, Interacting Particle Systems, (New- York: Springer). 

[19] T. M. Liggett, 1999, Stochastic Models of Interacting Systems: Contact, Voter and Exclusion 

Processes, (New- York: Springer). 
[20] C. Kipnis, S. Olla, and S. R. S. Varadhan, 1989, Hydrodynamics and large deviation for simple 

exclusion processes, Comm. Pure Appl. Math. 42, 115. 
[21] P. L. Krapivsky, S. Rcdncr and E. Ben-Nairn, 2010, A Kinetic View of Statistical Physics 

(Cambridge: Cambridge University Press). 
[22] A. J. Macfarlane, 1989, On q-analogues of the quantum harmonic oscillator and the quantum 

group SU(2) q J. Phys. A:Math. Gen 22, 4581. 
[23] K. Mallick, 2011, Some exact results for the exclusion process, Proceedings of Statphys 24, J. 

Stat. Mech.: Theor. Exp. P01024. 
[24] S. Prolhac, M. R. Evans, K. Mallick, 2009, The matrix product solution of the multispecies 

partially asymmetric exclusion process, J. Phys. A 42, 165004. 
[25] J. H. H. Perk and C. L. Schultz, 1981, New families of commuting transfer matrices in q-state 

vertex models, Phys. Lett. 84A, 407. 
[26] N. Rajcwsky, L. Santcn, A. Schadschneider and M. Schreckenberg, 1998, The asymmetric 

exclusion process: Comparison of update procedures, J. Stat. Phys. 92 151. 
[27] G. Sandow and G. M. Schiitz, 1994, On [7q[St/(2)]-Symmetric Driven Diffusion, Europhysics. 

Lett. 26, 7 

[28] B. Schmittmann and R. K. P. Zia, 1995, Statistical mechanics of driven diffusive systems, in 
Phase Transitions and Critical Phenomena vol 17., C. Domb and J. L. Lebowitz Ed., (San 
Diego, Academic Press). 

[29] C. L. Schultz, 1983, Eigenvectors of the Multi-Component Generalization of the six-vertex model, 
Physica 122A, 71. 

[30] G. M. Schiitz, 1997, Duality relations for asymmetric exclusion processes, J. Stat. Phys. 86, 
1265. 

[31] G. M. Schiitz, 2001, Exactly Solvable Models for Many-Body Systems Far from Equilibrium 

in Phase Transitions and Critical Phenomena vol 19., C. Domb and J. L. Lebowitz Ed., 

(Academic Press, San Diego). 
[32] E. R. Speer, 1993, The two species totally asymmetric exclusion process, in Micro, Meso and 

Macroscopic approaches in Physics, M. Fanncs C. Maes and A. Verbeure Ed. NATO Workshop 

'On three levels', Leuven, July 1993. 
[33] F. Spitzer, 1970, Interaction of Markov Processes, Adv. in Math. 5, 246. 

[34] H. Spohn, 1991, Large scale dynamics of interacting particles, (Springer- Verlag, New- York). 
[35] S. R. S. Varadhan, 1996, The complex story of simple exclusion, in Ito stochastic calculus and 

probability theory, 385 (Springer, Tokyo). 
[36] B. Wehefritz-Kaufmann, 2010, Dynamical critical exponent for Two-Species Totally Asymmetric 

Diffusion on a Ring, Proceedings of the XVIIIth International Colloquium on Integrable 

Systems and Quantum Symmetries SIGMA 6, 039 



